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Abstract 

We study the conservative and deterministic dynamics of two nonlinearly interacting particles evolv- 
ing in a one-dimensional spatially periodic washboard potential. A weak tilt of the washboard 
potential is applied biasing one direction for particle transport. However, the tilt vanishes asymp- 
totically in the direction of bias. Moreover, the total energy content is not enough for both particles 
to be able to escape simultaneously from an initial potential well; to achieve transport the coupled 
particles need to interact cooperatively. For low coupling strength the two particles remain trapped 
inside the starting potential well permanently. For increased coupling strength there exists a regime 
in which one of the particles transfers the majority of its energy to the other one, as a consequence 
of which the latter escapes from the potential well and the bond between them breaks. Finally, for 
suitably large couplings, coordinated energy exchange between the particles allows them to achieve 
escapes — one particle followed by the other — from consecutive potential wells resulting in directed 
collective motion. The key mechanism of transport rectification is based on the asymptotically van- 
ishing tilt causing a symmetry breaking of the non-chaotic fraction of the dynamics in the mixed 
phase space. That is, after a chaotic transient, only at one of the boundaries of the chaotic layer 
do resonance islands appear. The settling of trajectories in the ballistic channels associated with 
transporting islands provides long-range directed transport dynamics of the escaping dimer. 

PACS numbers: 05.60.Cd, 05.45.Ac, 05.60.-k, 05.45.Pq 



I. INTRODUCTION 



The study of transport phenomena has attracted considerable interest over the years due 
to its relevance in many physical situations. The latter are often described on the basis 
of one-dimensional particle motion in a tilted spatially periodic potential [lj-[13j. Corre- 
sponding experimental realisations include Josephson junctions [H], charge density waves 
[T5] . superionic conductors [16j, rotation of dipoles in external fields [TT] . phase-locked loops 
[T8j and diffusion of dimers on surfaces [H El [H HH1 [201 12U I221 [23] to name but a few. In 
many of these aforementioned situations the particles, in addition to their motion in the 
periodic potential, interact, which may lead to cooperative effects not found in situations of 
individual particle motion [24|-[27j. 

The objective of the current work is to investigate the conditions under which it is possible 
to generate a directed flow along with collective motion in a system of coupled particles. To 
be precise, we study the transport of a dimer evolving in a washboard potential experiencing 
a weak tilt force. The nonlinear bond dynamics between the two monomers, constituting 
the dimer, is modelled by a Morse potential allowing for bond breaking, i.e. fragmentation. 
We focus our interest on the chaos-promoted detrapping mechanism for dimers that initially 
reside in one well of the washboard potential. Provided that such a detrapping transition 
happens the question then is under which circumstances subsequent directed long-range 
particle transport is achievable. Since the total system energy is too low for both monomers 
to be able to escape from the potential well simultaneously, we explore whether coopera- 
tive energy redistribution is possible allowing at least one of the monomers to escape and 
subsequently display directed motion. We also elucidate the possible scenario in which the 
energy exchange between the monomers proceeds in such a well-coordinated manner that 
the monomers move separately from one well into the next, one following the other, resulting 
in directed motion of the dimer. 

The paper is organised as follows: In the next section the model of the dimer system is 
introduced, followed by the formulation of the escape problem together with a brief discussion 



of the related phase space structure. In Section III the particle current is studied and the 



occurrence of different transport scenarios is described. Afterwards in Section [TV] we relate 
the phase space dynamics to the regime of high particle current. In particular chaotic 
invariant sets, their connection with singularities of the escape time function, and their 
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relevance for the escape process are considered. In Section [V] we present an alternative 
description of the escape problem as the motion of a single particle in a two-dimensional 
potential landscape. Finally we summarise and discuss our results. 



II. THE MODEL OF THE DIMER SYSTEM 



We study the dimer dynamics with a Hamiltonian of the following form 



2 r 2 
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H int {q u q 2 ), (1) 



n=l L 

wherein p n and q n , n = 1, 2, denote the canonically conjugate momenta and positions of the 
two coupled particles of unit mass evolving in the periodic, spatially-symmetric washboard 
potential of unit period given by 

U M = U 0ig+ l) = -"^l. (2) 

The external field 

Ui{q) = ~F{q - log[cosh(g - q )\) (3) 

exerts a tilt on the washboard potential. The potential is sketched in Fig. 1 for tilt strength 
F = 0.01. The value of the parameter q = 10 in the second term on the right hand side of 
Eq. ([3]) is chosen such that the tilt rapidly diminishes when the coordinate q exceeds go and 
eventually upon further growth of q the bias vanishes. On the other hand as long as q -C qo 
the tilt adopts the value IF. Therefore particles that manage to escape from a potential 
well into the asymptotic region q < q — > 00 experience only a finite acceleration period 
at the end of which any forward motion must proceed unbiased. The question then arises 
whether escaping particles carry on moving forward even in the range where the bias is no 
longer present. 

The interaction part of the Hamiltonian Hi nt given by 

H^t = y(l - exp[-ct(g 2 -qi~ to)]) 2 , (4) 

is responsible for the coupling between the monomers which results from a Morse interaction 
potential of depth D, where a is the range parameter and the parameter l Q denotes the 
equilibrium distance between the monomers. Throughout the paper we chose Zo — 0.5, i.e. 
the equilibrium distance amounts to half the length of one period of the washboard potential. 



The equations of motion are 

q\ = - sin(27rgi) + F(l - tanh(gi - q )) 

+ aD(l - exp[-a(g 2 - qi - k)]) exp[-a(g 2 - qi - k)] (5) 

<?2 = - sin(27rg 2 ) + F(l - tanh(g 2 - qo)) 

- aD(l - exp[-a(g 2 - q 1 - l )]) exp[-a(g 2 - q 1 - l )] . (6) 

The interaction strength between the two monomers is effectively determined by the product 
aD. For aD = the system decouples into two integrable subsystems and the dynamics 
is characterised by individual regular monomer motions in the washboard potential. For 
nonzero aD the dynamics is no longer integrable. To prevent unphysical events in which the 
left monomer overtakes the right one, a sufficiently strong coupling between them is required. 
The choice a = 3 and D e [0.5, 3] ensures that. On the other hand the effective coupling 
strength, aD e [1.5, 9], is then too large by far to treat the coupling using a perturbational 
approach. We therefore resort to a numerical analysis of the coupled monomer dynamics. 

Let us briefly discuss the phase space structure corresponding to the dynamics in the 
tilted washboard potential. In the range — oo < q < 10 the tilt force, —dUi/dq, is effectively 
of strength 2F. For uncoupled monomers (aD = 0) there exist saddles at q k s = 0.5 + k — 
arcsin(2F)/(27r) and centers at q k c = k — arcsin(2F)/(27r) for integer values k. For very 
small tilt strength F < 0.01 the barrier height of the washboard potential, given by the 
difference between the energy of the saddle and the center, is virtually equivalent to those 
of the corresponding unbiased system with F — 0, i.e. E b ~ 1/ir. 

III. PARTICLE CURRENT 

In this section we consider the emergence of a particle current. The initial positions of the 
monomers are taken as —qi(0) = 92(0) = 0.25, so that the dimer is contained in one of the 
wells of the washboard potential and, for the weak tilt strength F = 0.01 used throughout the 
paper, is initially situated in virtually the lowest energy dimer configuration compatible with 
the bond length l . The dimer initially has potential energy E pot = J2 n =i 2t^o(?n) + Ui(q n )] = 
0.3221 which is of the order of the barrier energy E b ~ l/ir ~ 0.3183 of the washboard 
potential. Note that since the bond between the monomers is initially undistorted, the 
contribution from the Morse potential energy H int to the system's initial potential energy 

4 



is zero. The initial kinetic energy of the dimer is taken as Eun — 0.1234. Crucially the 
total energy E total = 0.4455 < 2E b is not sufficient that the two monomers can escape 
simultaneously from a well of the washboard potential. In order for directed motion of 
the dimer to occur at all, cooperation between the monomers, in the form of appropriately 
coordinated energy exchanges, is required. 

Particle transport is assessed quantitatively by the mean momentum, viz. the current, 
which is defined as the time average of the ensemble averaged momentum, i.e. 

P =^ r 'dt> \ Pl {f) + P2 (t')) , (?) 

J s Jo 

with simulation time T s and with the ensemble average given by 

1 - 

n=l 

Here N denotes the number of particles constituting the ensemble. For the computation 
of the ensemble average, trajectories belonging to iV = 2 x 10 5 values of the pair of initial 
momenta (pi(0),p2(0)) are taken. These initial values are uniformly distributed on an iso- 
energetic ring in the p\ — p 2 — plane such that the relation 

2E kin =p\+pl (9) 

is fulfilled. Notice the symmetry pi <-> — pi and i = 1,2. Hence there is no bias contained in 
the ensemble of initial conditions. The simulation time interval is T s = 10 5 being equivalent 
to almost 4 x 10 4 the period duration for harmonic oscillations near the bottom of a potential 
well. 

In what follows we vary the depth of the Morse potential, D, playing, for fixed a = 3, the 
role of the coupling parameter. The dependence of the current, defined in Eqs. ^ and 
on the value of D is shown in Fig. [2j For values D < 3.9 the current exhibits variations and 
even vanishes for D = 1.1. Interestingly the current effectively grows for d > 1.2. Finally the 
current rises rapidly and monotonically in the range 3.9 < D < 4.5 and effectively saturates 
at a high level for D > 4.5. We emphasise that the amplitude of the tilt force F = 0.01 is 
too small to alter the washboard potential significantly compared to the case without tilt. 
In fact the influence of the tilt force is sufficiently small that, for example, the potential 
barrier immediately to the right of the initial well is lowered by a mere 7.5% and hence the 
induced bias is very weak. 
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Figure 2: Current as a function of D in dimensionless units. The remaining parameter values are 
given by: a = 3, F = 0.01 and qo = 10. 

In the following we illustrate the complex solution behaviour of system (J5]) , ([6]) and the 
implications for the contribution to the net current. Fig. [3] shows the temporal evolution of 
the coordinates qi(t) and g 2 (t) for four different values of D but for the same initial condition 
leading to various types of solutions. For the low value D = 0.5 the coordinates perform 
small-amplitude oscillations around their respective starting value (see the upper left panel 
of Fig. [3J. Thus the monomers remain trapped in the potential well and the contribution 
to the net current is zero. In contrast for D = 1 we observe that after a finite period of 
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Figure 3: Four qualitatively different scenarios, illustrated by the time evolution of individual 
coordinates with initial conditions pi(0) = —0.2100 and pa(0) = 0.4502 for the values of D indicated 
in the plots. The remaining parameter values are given by: a = 3, F = 0.01 and qo = 10. The 
insets show details of the evolution. 

chaotic but bounded dimer motion the bond between them breaks. As a result the right 
monomer (with index n = 2) is released and due to the (still acting) tilt force accelerated 
into the region of higher momenta while its left counterpart (with index n = 1) becomes 
again trapped in a potential well. The subsequent regular dynamics is characterised by 
different motions of the monomers, namely that of the right monomer moving rightwards 
(rotations) in the asymptotic region and the left monomer performing bounded oscillations 
in a potential well (librations). In this case, the directed motion of the right monomer gives 
a contribution to the net current. We stress that, after such fragmentation, reformation of a 
bound state dimer from the two isolated monomers is excluded. Notice that the possibility 
of bond breaking allows for transient chaos [28]. [29]. [30]. 
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Increasing the coupling parameter further to D = 2 the chaotic coupled monomer dy- 
namics involves irregular phases where the motion changes from forward to backward and 
vice versa in the manner of Levy flights |31| in the whole simulation interval. Nevertheless 
the net motion proceeds to the right. Notice that the bond remains unbroken, and hence 
the dimer intact, in this case. 

Interestingly for a high value D = 3 the dimer manages quickly to escape from the po- 
tential well. (We remark that for some other initial conditions we observed first a longer 
transient of still bounded but chaotic motion before the escape eventually took place.) Fur- 
thermore, as the inset reveals, the two monomers perform out-of-phase motion, viz. the 
length of the bond between them alternately (slightly) decreases and increases. This is as- 
sociated with such well-coordinated energy exchanges between the monomers that firstly 
the right-hand monomer overcomes the potential barrier and reaches the adjacent well on 
the right, subsequently the left monomer follows, and so on. We underline that prior to its 
arrival in the asymptotic region the trajectory passes through a chaotic transient to adopt 
regular dynamics in the asymptotic region. It is the asymptotic vanishing of the tilt force 
that makes transient chaos possible. Clearly the directed dimer motion contributes with 
significant weight to the net current. In particular, the dimer moves with higher velocity for 
D = 3 than the escaped monomer does for D = 1. 

To summarise briefly: we distinguish between four qualitatively different transport sce- 
narios: 

(i) The dimer remains trapped inside the starting potential well and hence, there results no 
contribution to the net current. 

(ii) The dimer escapes from the starting potential and undergoes subsequently diffusive mo- 
tion with no substantial contribution to the net current. 

(iii) Directed energy transfer from the left monomer to the right one leads not only to frag- 
mentation but also to such a high energy gain of the right monomer that it can undergo 
directed motion to the right. This individual directed motion yields a considerable contri- 
bution to the net current. 

(iv) Appropriately coordinated energy redistribution between the monomers leads to di- 
rected collective motion such that first the right monomer performs a transition from one 
potential well into the next one to the right and afterwards the left monomer follows. Notice 
that this corresponds to repeated detrapping-trapping transitions. This scenario is optimal 
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in the sense that it enables both monomers to perform consecutive step-wise escapes from 
the potential wells. Since the total amount of energy does not suffice for a simultaneous 
escape of both monomers they must necessarily share their energy cooperatively in order 
to achieve escape at all. This almost-periodic energy exchange between the monomers cor- 
responds in phase space to motion near a stable period-one fixed point (see further in |V|). 
Notably the directed chaotic motion persists even in the asymptotic region where there is 
no bias anymore. Furthermore, the resulting velocity is higher than in case (iii) and so is 
the contribution to the net current. We mention that, apart from the ideal situation of on- 
going directed motion, for other initial conditions the dimer performs directed long-distance 
motion in a restricted time interval at the end of which non-directed diffusive motion as in 
case (ii) follows. 

We emphasise that the scenarios shown in Fig. [3] are not necessarily representative of 
the dynamics of all initial conditions at the respective values of the coupling strength D; 
at D = 3, for example, scenarios (iii) and (iv) both occur (for different initial conditions). 
Concerning the current we therefore remark that, at each fixed value of the coupling strength 
.D, each of the transport scenarios (i)-(iv) will, if present, contribute with different weight to 
the ensemble average for the current, which results in the complex behaviour seen in Fig. [2} 

IV. PHASE SPACE DYNAMICS 

In order to discuss the corresponding dynamics taking place on the three-dimensional 
energy hypersurface in the four-dimensional phase space we introduce the following Poincare 
surface of section (PSS) 

S = {p 2 ,g 2 |gi = 0,p 1 >0}. (10) 

In Fig. 4 (a) and (c) we depict the PSS for D = 0.5 and D = 3, respectively using an ensemble 
of 10 4 initial conditions fulfilling the relation ^ with E kin = 0.1234. The corresponding 
right panel presents the escape time function, defined as the time it takes the right monomer 
to reach the position q 2 = 10, as a function of the angle 

$ = tan- 1 (p 2 (0)M(0)). (11) 

For D = 0.5 there are two wide regions on the $— axis for which no escape happens at 
all. This is the case when the monomers start with initial momenta of approximately equal 
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Figure 4: PSS represented in the p 2 — plane (a) and (c) and escape times versus the angle 
= tan -1 (j>2(0)/pi(0)) shown in (b) and (d). Upper (resp. lower) row: D = 0.5 (resp. D = 3). 
The remaining parameter values are given by: a = 3, F = 0.01 and qo = 10. An ensemble of 10 4 
initial conditions fulfilling relation ^ with Eki n = 0.1234 is used. 

magnitude but of opposite sign, i.e. around $ ~ 3n/A and $ ~ 7n/A. In both cases the 
resulting regular motion is associated with the stable island centered at (p 2 , q 2 ) — (—0.2, 0.4) 
in the corresponding PSS (displayed in the inset of Fig. [31(a)). The physical reason for the 
appearance of regular trapped motion in these cases is the fact that the initial velocity of the 
center of mass of the dimer is virtually zero — a situation that remains virtually unchanged 
due to the symmetry of the washboard potential, the weakness of the tilt, and the small 
interaction term. An examination of the escape time function at various scales reveals that, 
except for the two windows of no-escape, the escape time depends sensitively on changes of 
the initial momenta. The trajectories attributed to escaping monomers are contained in the 
extended chaotic sea which densely fills the majority of the energetically-accessible regions 
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on the PSS (except for an infinite set of smaller islands of stability not recognisable on the 
scale of the PSS). A crescent-moon-shaped region within the chaotic sea remains empty 
on the PSS because it is not energetically accessible. In more detail, there exist chaotic 
invariant sets consisting of homoclinic and heteroclinic tangles which induce a fractal set of 
singularities into the escape time function (32]- (36] . The singularities arise at those points 
where the stable manifolds of unstable periodic orbits intersect the set of initial data with 
the effect that the corresponding trajectories become trapped for arbitrarily long times in a 
chaotic invariant set. It is therefore impossible to fully resolve the behaviour of the escape 
time function whose singularities form a fractal set with measure zero. 

Interestingly, for increased coupling strength D = 3 the interaction between the monomers 
is strong enough that fully developed chaos results and the windows of no-escape obtained 
in the previous case of D = 0.5 vanish. In comparison, the escape times are mostly shorter 
by far for D = 3 than for the preceding case D = 0.5. The associated PSS elucidates these 



differences in the escape process. (We remark that according to the condition in (10) only 
those trajectories for which the left monomer is still in the starting potential well contribute 
to the PSS.) Comparing the cases D = 0.5 and D = 3 one infers that in the former case the 
chaotic sea engulfs far more area (extending along the q 2 — axis over the range of the starting 
potential well together with its neighbour to the right) despite the existence of the (small) 
stable island (inset in Fig. [4](a)). Moreover the fact that the PSS is more densely populated 
for D = 0.5 than in the case D = 3 indicates that trajectories spend longer times in the 
potential well(s) for D = 0.5 before they manage to escape. Furthermore, for D = 3 some 
trajectories follow directly the unstable manifold associated with a chaotic saddle appearing 
as a winding curve emanating from the region around p 2 = 0,q 2 — 0.5 (inset in Fig. [3](c)). 
Further details are given in Section |V} This provides a mechanism for fast escape into the 
range of large coordinates which happens particularly for initial values lying in the range 
< $ < 1.5. Nonetheless there remains a large portion of trajectories that dwell for some 
time in the starting region before they escape in the direction of the asymptotic region. The 
dwell time depends sensitively on the initial conditions. On the other hand, escape does not 
necessarily imply sustained directed transport. 

For further characterisation of the escape dynamics PSS are presented in the pi—q\ — plane 
using the intersection condition q 2 = 20, p 2 > yielding a snap shot of the ensemble dynamics 
in the asymptotic region. Note that at such a large distance from the starting point located 
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Figure 5: Sections taken when the coordinate of the right monomer reaches q 2 = 20 with 10 4 initial 
conditions from an ensemble satisfying the relation ^ with E^ in = 0.1234. (a) (resp. (b)): D = 0.5 
(resp. D = 3). The remaining parameter values are given by: a = 3, F = 0.01 and qo = 10. 

at 52(0) = 0.25, the right monomer experiences an unbiased potential. The ensemble of 
initial conditions is the same as the one used for Fig. [4| The PSS for D = 0.5 is depicted 
in Fig. [5](a) and shows that advancing right monomers leave the overwhelming majority of 
their left counterparts behind distributed over various potential valleys where they perform 
trapped motion. Clearly for the fairly low interaction potential depth D = 0.5 the bond 
between the monomers easily breaks. In contrast, for the comparatively large interaction 
potential depth, D = 3, the bond between the monomers remains intact for the entire time, 
as seen in the right panel of Fig. [5j Hence, left monomers travel the full distance in unison 
with their right counterparts. 

It is illustrative to consider the distribution of the momenta of the right monomers at 
the moment when they reach the position q 2 = 20. For D = 0.5 the momenta are narrowly 
distributed around the peak value p 2 — 0.92 (not shown). In this case the particle transport 
is dominated by directed motion of right monomers after fragmentation (cf. scenario (iii) 
above). For D = 3 there results a broad momentum distribution in an interval matching 
that covered by the p\— values in the right panel of Fig. [5j This indicates that the dynamics 
in the asymptotic region involves not only directed motion but also itinerant motion as 
described above by the diffusive-like scenario (ii). Nevertheless the distribution of the p 2 
values attains a maximum at p 2 — 0.81, viz. for momenta for which directed motion to the 
right proceeds. 
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V. MOTION OF A PARTICLE IN AN EFFECTIVE TWO-DIMENSIONAL PO- 
TENTIAL 



To gain further insight into the coupled monomer dynamics it is useful to perform the 
following canonical change of variables induced by the generating function: S = \{qi + 
Q2)Px + \{<l2 — Qi)Py relating the old and new variables as follows 

Pl = \(Px + Py), P2 = \{P x -Py), (12) 

Qx = qtei + ftO, Q y = 2^2 - qi) ■ (13) 

The coordinate Q x determines the position of the center of mass (CM) of the dimer, ac- 
counting for translational motion. Vibrations (V) of the dimer are described by Q y . The 
Hamiltonian expressed in the new variables becomes 

H = -(Pi + P„ 2 ) - - cos(2ttQ x ) cos(2nQ y ) 

4 y 7T 

+ ^(l- exp [-a(2g y -/ )]) 2 

- F(2Q X - log[cosh(g a -Qy- q )\ - log[cosh(Q a; + Q y - q Q )\) (14) 
= \{F* + I*) + U{Q m ,Q y ). (15) 

The corresponding equations of motion, describing the effective motion of a particle in a 
two-dimensional potential landscape U(Q x ,Q y ), are given by 

Qx = —2sm(2TrQ x )cos(2iiQy) 

+ F[2 - tanh(Q x - Q y - q ) - tanh(Q x + Q y - q )] 

Q y = —2 cos(2ttQ x ) sm(2iiQy) 

- 2aD{\ - exp[-a(2Q y - l )}) exp[-a(2Q y - l )} 

+ F[tanh(Q x - Q y - q ) - tanh(Q a . + Q y - q )] . (16) 

For Q X iQy <C go the impact of the external tilt force matters only in the first equation 
whereas the Morse coupling enters only in the second equation. The interaction between the 
Q x (CM) and Q y (V) degree of freedom (d.o.f.) results from parametric modulations of the 
respective washboard potential force term. The effective potential U(Q X , Q y ) is displayed 
in Fig. [6] for the interaction potential depths D = 0.5 and D = 3 respectively. The super- 
imposed trajectory, corresponding to the dynamics shown in Fig. |3](a) and (c) respectively, 
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Figure 6: (Colour online) Evolution of the trajectory in the two-dimensional potential energy 
landscape U(Q x ,Q y ). The parameter values are Iq = 0.5, a = 3. Left panels (a,c,e): trapped 
particle for D = 0.5 and right panels (b,d,f): moving particle for D = 3. For clarity, the middle 
(c,d) and bottom (e,f) rows show profile views (Q X ,U) and plan views {Q x ,Qy), respectively. The 
steepness of the potential surface in the case D = 3 (b,d,f) necessitates plotting the potential surface 
for a slightly smaller Q y range. The left and right panels are shown to the same scale in each case. 

starts close to the potential minimum for D = 0.5 situated at (Q X1 Q y ) = (0.014, 0.173). The 
corresponding state of lowest energy is denoted by E g . There exists a nearby saddle, afore- 



mentioned in Section IV, which for D = 0.5 is located at (Q x ,Q y ) = (0.241,0.244) having 
energy U s . In order to advance towards higher Q x — values in the two-dimensional potential 
landscape the particle needs to overcome a potential barrier the height of which is deter- 
mined by AU — U s — U g . Apparently for D = 0.5 the trajectory remains trapped inside the 
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Figure 7: (Colour online) Difference between the energy of the potential minimum and the nearby 
saddle as a function of D. The remaining parameter values are Iq = 0.5, qo = 10, F = 0.01, and 
a = 3. 

potential well (cf. Fig. |3](a)). This is mainly connected with the relatively large-amplitude 
excursions in the Q y — direction pointing to rather pronounced bond stretching/compression. 
In fact, the major part of the total energy is contained in the Morse interaction term, viz. 
in the V-d.o.f., amounting to 70%. Thus there remains little energy that can flow into the 
CM-d.o.f., hampering the translational motion necessary to overcome the potential barrier. 
In contrast for D — 3, when the bond between the monomers is more rigid by far than 
before, the Morse bond energy constitutes only a small amount of the total energy. As a 
consequence the CM-d.o.f. possesses enough energy that the trajectory easily overcomes 
the potential barrier and passes from one well to a neighbouring one (right panel in Fig. [6]). 
Moreover, along the line Q y = 0.25, that is —q\ = q 2 = 0.25, there is no gradient of the 
potential in the Q x — direction (CM motion direction). Therefore a strong enough interac- 
tion strength D is advantageous for transport because it confines the motion of the dimer 
along a narrow strip centered along the line Q y = 0.25. At the same time the height of 
the energy barrier AU decreases with increasing interaction potential depth as illustrated 
in Fig. [7j Conclusively, motion of the mean coordinate Q x from one potential well into the 
neighbouring one is readily accomplished for large values of D which is reflected in a high 
current (see Fig. [2] above). 
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Figure 8: (Colour online) PSS for the dynamics corresponding to a symmetrical uniform distribution 
of initial conditions satisfying relation ^ in the potential well at Q x = 0, showing (a) Q x < 10 
and (b) the asymptotic regime Q x > 10 (shown mod (1)), for coupling parameter D = 3. The 
remaining parameter values are a = 3, F = 0.01 and qo = 10. The panel (c) shows detail of the 
mixed phase space at the upper boundary of the chaotic layer in (b). 

Finally, we relate the escape process and the emergence of directed chaotic motion to the 



phase space structure of the transformed system with Hamiltonian given in Eq. (15). To 
this end we use the following PSS 



Z = {P x ,Q x \Q y = l/4,P y >0} . 



(17) 



In Fig. [8] we plot for a strong particle coupling D = 3 the PSS corresponding to the escape 
dynamics when Q x < 10 and the dynamics in the asymptotic region, i.e. Q x > 10 (for 
those trajectories which reach it), in (a) and (b) respectively, being characterised by chaotic 
sets. (In Fig. |8](b) the coordinate Q x is presented mod(l).) In Fig. |8](a) chaotic saddles, 
formed by the intersecting stable and unstable manifolds of unstable periodic points, govern 
the dynamics. The majority of escaping trajectories follows the unstable manifold of the 
saddle point located at (Q x ,Q y ) = (0.188,0.243). On the other hand there are trajectories 
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that remain in the starting region or spend at least some time there before escape as a 
consequence of the presence of chaotic saddles , [3H] - Furthermore, on approaching the 
asymptotic region, where the tilt of the washboard potential vanishes, some of the previously- 
escaping dimers become trapped in wells of the washboard potential again. 

From Fig. |8](b), depicting the PSS in the asymptotic region, we conclude that the bulk of 
the layer on the PSS is covered by a chaotic set. Within the chaotic set trajectories move in a 
diffusive way with changes of the direction of motion not contributing to transport. Notably, 
at the upper boundary of the layer, shown in Fig. j8](c), islands of regular motion arise 
from those trajectories that have settled on regular dynamics after their passage through a 
chaotic transient. Most importantly, these islands possess non-zero winding numbers and 
thus act as ballistic channels [39] providing directed transport to the right. In particular, the 
dynamics within the island structure centered at the stable period-one fixed point (P x , Q x ) = 
(0.703,0.180) reflects the almost-synchronous monomer motion described in scenario (iv) in 



Section III In more detail, motion near the fixed point corresponds to almost-periodic energy 
exchange between the monomers which is connected with only minor bond deformations 
which corroborates the findings reported above for the optimal transport scenario. 



VI. SUMMARY 



We have analysed the Hamiltonian dynamics of two nonlinearly coupled particles evolving 
in a washboard potential. Notably the total energy does not suffice to enable simultaneous 
escape of the two particles, initially trapped in a well of the washboard potential. Due 
to appropriate energy redistribution, at least one of the particles can achieve escape. (See 
also [10].) Ideally the two particles share energy almost periodically in such a way that 
consecutive detrapping-trapping transitions take place during which the particles escape 
one after another from one well into an adjacent one. It has been demonstrated that a 
weak tilt force, vanishing asymptotically in the direction of the bias, is sufficient to instigate 
directed motion of the escaping particles. Transport is accomplished for those trajectories 
which follow a chaotic transient, associated with the dynamics of chaotic saddles, settling 
afterwards on regular motion. The key mechanism of current rectification is based on the 
asymptotically vanishing tilt causing a symmetry breaking of the non-chaotic fraction of the 
dynamics where only at the upper boundary of the chaotic layer resonance islands appear. 
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The latter is supported by transporting island structures in the mixed phase space which 
serve for long-range directed transport. 

Finally, we mention that it is certainly interesting to extend the study of directed motion 
to systems involving many more degrees of freedom than in the current dimer case where the 
dynamics within transporting islands can be investigated utilising two-dimensional Poincare 
surface of sections. In particular, it needs to be examined what structures in higher dimen- 
sional phase spaces play the role of possible ballistic channels providing directed collective 
transport. 
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